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Abstract: In this paper we show that the sensitivity of absorption 
imaging of ultracold atoms can be significantly improved by imaging in 
a standing-wave configuration. We present simulations of single-atom 
absorption imaging both for a travelling-wave and a standing-wave imaging 
setup, based on a scattering approach to calculate the optical density of a 
single atom. We find that the optical density of a single atom is determined 
only by the numerical aperture of the imaging system. We determine 
optimum imaging parameters, taking all relevant sources of noise into 
account. For reflective imaging we find an improvement of 1.7 in the 
maximum signal-to-noise ratio can be achieved. This is particularly useful 
for imaging in the vicinity of an atom chip, where a reflective surface is 
naturally present. 
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1. Introduction 

The detection of single atoms with high fidelity is an important requirement for many modern 
atomic physics experiments, and in particular for quantum information science. Single particle 
imaging has been achieved using fluorescence imaging of ions [ 1 1 and neutral atoms |l2}{6l, 
but is difficult in atom-chip experiments due to a scattering background from the chip H. 



A different approach is to enhance the atom-photon interaction using resonant cavities |l7H9l. 
However, free-space absorption imaging of single particles lfT0HT2l is still a daunting task. 

In this paper we propose absorption imaging in a standing wave as a simple and effective way 
to improve the signal-to-noise ratio (SNR) of atom detection. This is particularly useful in the 
context of atom-chip experiments 1131 where a standing wave naturally forms if the imaging is 
done perpendicular to the chip surface. We will take the setup currently in use in our group as 
a reference [14-16| for realistic experimental parameters. 

In the first section we will calculate the absorption signal of a single atom, both for traditional 
one-pass imaging and for reflection imaging as employed in our experiment. In the second half 
of the paper we present numerical simulations of the expected absorption images and signal- 
to-noise ratio (SNR) under realistic conditions, taking into account all common sources of 
noise, for in-trap atoms as well as for untrapped atoms. To do so we use a hybrid approach, 
first calculating the movement of the atom in the imaging beam based on scattering individual 
photons, and subsequently investigating the absorption signal due to this atom, now taking the 
imaging beam as a classical wave. We show that using reflection imaging drastically increases 
the expected SNR and thus the single-atom detection accuracy. 

2. Single-Atom Absorption Imaging 

Our starting point is the absorption signal of a single atom at rest, illuminated by a resonant 
laser beam with homogeneous intensity distribution (see Fig.[T]i. As is shown in the appendix 
the intensity in the imaging plane is then given by 

^^^[Rewp))-„|.(P)P] a, 

where /abs(p) is the absorbed intensity, i.e. the missing intensity in the image plane with image- 
plane coordinates p, 1^ is the incoming intensity, a is the absorption cross section of the atom, 
s the saturation parameter in the object plane, p{p) the point-spread function of the imaging 

system (see below), a — [J \p{p)\^ d^p] the effective area over which the signal is distributed 
in the detection plane and x is the fraction of the scattered light collected by the imaging system. 

In an experiment the above quantities a and /abs (0) are not directly accessible due to the finite 
size of any detector and the finite time resolution. Instead one would observe the quantities 
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where x is an index over camera pixels, A is the size of a pixel and T is the exposure time, i.e. the 
duration of the imaging pulse. Both a and s now depend on time, as the atom acquires a finite 
velocity during the imaging pulse, leading to Doppler shifts. In addition the local saturation 
parameter in the object plane varies in the case of standing wave imaging, as will be discussed 
below. We can then define the apparent column density of atoms per pixel in the image plane 
as 

n{x) ^ — — , . (3) 

where sq is the saturation parameter for the intensity of the incoming beam and Gq is the ab- 
sorption cross section of an atom at rest for the incoming intensity. The apparent atom number 
per pixel would be given by Napp{x) =An{x). 



As the atoms will move out of the focal plane during the imaging pulse, simulating this signal 
will require a point-spread function that is a function not only of the image plane coordinates, 
but also takes defocussing into account. This PSF can be numerically determined as 
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p{p,f;t) = / rexp(//r^) x Jq {2Kpr/p^)dr (4) 



for any given defocus / = ^'Z.{1 — V 1 — NA^) where Z is a real-space coordinate in the imag- 
ing direction 1 17 1, NA is the numerical aperture of the imaging system and both p and Z might 
depend on time, po = A /NA is proportional to the resolution of the imaging system, determined 
by the wavelength and the numerical aperture. Aberrations in the imaging system can be taken 
into account in a similar manner. We assume these to be negligible near the optical axis. In the 
simulations the numerical point-spread function given above is used to calculate the intensity 
in the image plane. 

2.1. Analytical Estimates 

For a better intuition of the above results we consider some simple analytic estimates. Assuming 
a stationary atom in perfect focus of a two-lens imaging system with unit magnification (/-2 /- 
/) we can use a simple point-spread function of the form p{p) = l/(ppo)/i(27rp/po) where 
^(O) = tt/Po = Using x — (V8)NA^ and assuming a 2-level system where ao = 3X^/2n 
the maximum signal at the center of the image simplifies to 

hn{0) l+so \2 16 / 

which is thus determined purely by the numerical aperture of the imaging system and the sat- 
uration parameter sq. For a numerical aperture NA = 0.4 as is the case in our experiment we 
then find an absorption in the centre of the image of /abs(0)//in(0) — 0.23 in the low saturation 
limit 50 ^ 1 for a single atom. 

We can further consider the total apparent number of atoms one would extract from this 
signal as the peak signal times the area over which the signal is distributed. 



,l+50/abs(0) 
do /in(0) 



^app =a^^^^ = 1 - -NA2 (6) 



This shows that the expected number of atoms is reduced by the second, NA-dependent term. 
This term is due to light which is scattered by the atom, but scattered into the solid angle of the 
lens, thus reducing the apparent amount of light that is absorbed. For a numerical aperture of 
0.4 we therefore expect to find an apparent atom number A^app of 0.94 rather than one. 

The commonly used quantity optical density (OD) can thus be defined for a single particle 
by the point-spread function of the imaging system, i.e. OD = — ln(l — 4bs//i„). At the center 
of the image (for vanishing pixel size) we find for the peak optical density a value of OD(0) ~ 
3/2NA^ + ^/leNA"* in the low-saturation limit (where the sign change is due to a series expansion 
of the logarithm). For our numerical aperture this equals a peak optical density of 0.26. It is 
worth noting here that the atom number extracted in the usual way from optical density, in this 
case 



would yield 1.07 atoms rather than 0.94. The atom number is slightly overestimated, as 
Lambert-Beers law is valid only for a continuously absorbing medium, not for a single ab- 
sorber In the following we always use the spatially dependent equivalent of (|6]l rather than (|7]i 
for the atom number 



2.2. Reflection Imaging 
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Fig. 1. Sketch of the imaging system considered for this paper. On the left the situation 
of a travelhng wave without reflection is depicted. On the right the situation of reflective 
imaging (with the atom chip on top) is shown. 

The above discussion is for an atom imaged by a travelling wave, as depicted in Fig. [ija). 
In reflection imaging, such as in Fig. [TJb) the situation is slightly different. In this case the 
probe light forms a standing wave at the position of the atoms. For Rb atoms in our magnetic 
lattice lfT4l4T6ll we expect a ground-state size of the atomic wavefunction of about 40 nm, about 
a factor of 10 smaller than the period of the standing wave. We therefore approximate the atom 
as a point-particle which can initially be positioned in an anti-node of the standing wave. In our 
atom chip experiment the trap position can be adjusted by means of an externally controlled 
magnetic field to achieve this. 

The intensity of the probe beam then varies sinusoidally along the probe direction and the 
maximum intensity in the anti-nodes is a factor four greater than the intensity of the incoming 
travelling wave. The coherent scattering amplitude £/ scales with the the local field amplitude 
(at low saturation). Thus the scattered wave is two times as strong as for a travelling wave, 
leading to an accordingly higher observed signal. More precisely, the saturation parameter will 
also change. 

It should be noted here that while in the limit of low saturation and short exposure times the 
observed signal will indeed be two times larger than for a travelling wave, for realistic imaging 
with finite exposure time the atom starts to probe the spatial varation of the light field, leading 
to an observed signal which depends on the exact imaging parameters. The expected signal can 
then only be predicted by simulations such as described below. 

A further important difference to travelling wave imaging is that the atom is not pushed out 
of the focus of the imaging system by the probe beam, but performs a random walk in all 
dimensions. This point wiU also be further discussed in the simulations below. 

Another factor of two could in principle be gained because the scattered wave is emitted 
towards and reflected by the mirror. This again doubles the amplitude of the scattered wave 
in the detection plane. A practical consideration here is that for single-atom imaging one will 
typically use high-NA imaging, with a Rayleigh length zr = X/ (ttNA^) < 10|J.m. For the 
typical situation in our experiment the atom is at about 10 \xm from the surface. This implies 
that the atom and its mirror image are not simultaneously in focus in the image plane. In practice 
this reflected wave hardly contributes to the optical density for our experimental parameters. If 
the imaging is done extremely close to a reflecting surface, or if the Rayleigh length of the 



imaging system is significantly larger this can however be a significant effect. 



3. Simulations 

After having calculated the absorption signal of a single atom at rest in the previous section, 
we would now like to determine the accuracy with which we can detect a single atom in a more 
realistic setup. To do so we have to take two important effects into account: recoil blurring and 
imaging noise. The first describes the process where the atom starts to move due to scattering 
photons, leading to a blurring of the image. The second describes the uncertainty with which 
one can determine the local intensity of the imaging beam, both for fundamental reasons, i.e. 
shot-noise, and technical reasons such as read-out noise of the camera. These two cases impose 
counteracting constraints: recoil blurring is minimal for short exposure times and weak probe 
intensities, shot-noise is reduced for large photon counts, i.e. long exposure times and high 
probe intensities. In addition, at high intensities the signal saturates. Therefore there exists an 
optimum signal-to-noise ratio (SNR) for a certain exposure time and probe intensity. We use 
a hybrid approach in treating these effects: recoil blurring is calculated based on scattering 
individual photons off the atom (see below), while the absorption signal is calculated based on 
a classical-wave approach as described above. 
In our simulations we distinguish four cases: 

a) free atom in travelling wave 

b) trapped atom in travelling wave 

c) free atom in standing wave 

d) trapped atom in standing wave 

where the cases a) and b) are for traditional one-pass imaging, and the cases c) and d) are for 
our reflective atom-chip setup (see Fig.[T]l. While case a) and c) assume a freely moving atom, 
case b) and d) assume the atom to be held in a magnetic trap. For this trap we assume the poten- 
tial calculated from the magnetic pattern of microtraps present on our atom chip 1 14| which is 
approximately equivalent to trap frequencies of (0)^,0)3,0-) = 27rx(38.1,36.5, 14.0)kHz and 
a trap depth of 8.8 G. In the simulations the calculated potential is used without harmonic ap- 
proximations. The atom is treated as free as soon as it acquires sufficient energy to escape the 
trap. We assume optical pumping to non-trapped states is negligible for an appropriate choice 
of polarisation for the imaging beam. 



case 


z(p.s) 


///sat 


SNRpx 


SNRcRB 


Napp 


PlA 


Pl.O 


a) 


12.4 
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1.73 
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0.74 


62% 
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d) 


42.5 


0.56 


2.63 


2.68 


0.83 


82% 


9.5% 



Table 1. Optimum exposure parameters, and resulting SNR for the four cases described in 
the text. SNRpx is the SNR obtainable from evaluating a single pixel, while SNRcrb is the 
SNR for an estimator achieving the Cramer-Rao bound (without taking fringe-removal into 
account). The normalization factor N^pp for case a) and b) is as expected from the analytic 
treatment of section [2T| The normalization factor in cases c) and d) can only be obtained 
from the simulations. Finally we give the probability of true positive measurements pn 



and false positive measurements due to the presence of zero atoms pi Q (see 3.2 1 



3.1. Recoil Blurring 

As the atom is illuminated by the imaging beam it scatters photons at a rate determined by the 
local intensity at the position of the atom as well as any detuning due to e.g. Doppler shifts as 
the atom acquires momentum with each of these scattering events. In the case of a travelling 
wave the photon scattering drives the atom strongly along the direction of the imaging beam, as 
the absorption of a photon always happens in this direction, performing a biased random walk. 
This quickly leads to the atom moving out of the focal plane of the imaging lens. The atom 
further performs an unbiased random walk in momentum space in the plane perpendicular to the 
imaging direction due to the spontaneous emission of a scattered photon, leading to a blurring 
of the image. In the case of a standing wave the atom can absorb a photon from both directions 
along the imaging axis, giving rise to an unbiased random walk in the imaging direction as well 
as the plane perpendicular to it. The atom therefore leaves the focus of the lens much more 
slowly than in the case of a travelling wave. 

We simulate the movement of the atoms by calculating stochastic trajectories where random 
scattering events occur at the local scattering rate, taking the effect of the trap into account as 
a classical potential. In the case of an untrapped travelling wave we recover the expected be- 
haviour, i.e. the root-mean-square (RMS) position along the imaging beam growing as and 
the RMS position in a plane perpendicular to the imaging direction growing as f^/^, as is ex- 
pected for uniform acceleration and a random walk in momentum space respectively. The effect 
of the trap is almost negligible in the case of a travelling wave, as the force of the light field is 
much greater than the restoring force of the trap potential. In the case of reflective imaging the 
behaviour proportional to is also expected in the probe direction. This is however modified 
by the periodic intensity of the standing wave. Since the net force on the atoms is much lower 
than in the case of a travelling wave the influence of the trapping field is also more important, 
particularly at low intensities. 

These calculated RMS trajectories are then used in the integrations of Eq. (|2]i to determine 
the signal after a given exposure time for a certain probe intensity. 

3.2. Detection Noise 

To determine the feasibility of absorption imaging for single atoms one needs to consider vari- 
ous sources of noise in addition to the atomic trajectories to determine the signal-to-noise ratio. 
In the following we assume near shot-noise limited imaging, with further contributions from 
camera read-out noise and dark counts, as well as a finite quantum efficiency of the imaging 
system. Equations (|2]i determine the number of photons per pixel in the signal and light fields, 
and can be trivially modified to take a finite quantum-efficiency of the imaging system into 
account. With these values we can calculate the expected signal n(x) as atom density per pixel 
(Eq.|3]l for any combination of exposure time and intensity. The variance of this signal is given 



where A'det and A^ief are the number of electrons per pixel in the signal and reference images 



of these quantities, mostly determined by photon shot noise (o^ A^) and camera readout noise. 

The signal to noise ratio per pixel is then given by «(^)/c7„(.i); one can further improve the 
signal to noise ratio by not only evaluating the central pixel, but using an optimal estimator 



by 




(8) 



respectively, also integrated over pixel size and exposure time. O'^ 



ilgj and ffj^f are the variances 



achieving the Cramer-Rao bound, such as 



N=^^ (9) 

where q{x) is the spatial mode function of the signal ifTSll . In our results we list the SNR both 
for evaluating only the central pixel and for an estimator achieving the Cramer-Rao bound. 

We define the accuracy of the measurement as the probability of finding true positives pi i, 
i.e. finding one atom if there is one atom present. Here we assume all measurement outcomes 
for which 0.5 < N < 1.5 to indicate the presence of exactly one atom. We have determined 
this value by simulations of 10*' individual absorption images for each case. The results are in 
excellent agreement with those expected for a normal distribution, for which the accuracy is 
given hy F = Erf{z/^) with z = ^/la. 

We can also determine the probability of finding false positives in the absence of any atoms 
pi,Q, i.e. finding one atom if zero atoms are present. False positives due to the presence of more 
than one atom can not easily be determined, as an accurate calculating of the absorption signal 
of two atoms is non-trivial. We expect these to be of similar magnitude as those due to the 
absence of any atoms. 

3.3. Simulation Results 



Trapped 




20 40 60 80 20 40 60 80 100 

Exposure Time {\is) 



Fig. 2. Simulated SNR as function of intensity of the incoming beam (in units of saturation 
intensity) and exposure time. The four subplots correspond to the four cases described in the 
text as indicated near the top and right axes, the color scale is normalized to the maximum 
SNR of case d). The black dot marks the position of the optimum. 



We simulate the expected absorption image and SNR assuming a numerical aperture of the 
objective lens of 0.4 (Edmund Optics NT47-727 with custom AR coating), a total magnification 
of the imaging system of 10, a pixel size of the camera of (13 x 13) \xm^ (Andor iKon-M 934), a 
quantum efficiency of 0.9 and a readout noise of 13 counts per pixel ifTSl . We consider resonant 
CJ+-transitions on the D2-line of Rubidium 87, where Co = 2.907 x 10^^ cm^ and the saturation 
intensity /sat — 1 .669mW/cm^ 1 19 1. We have confirmed that in the limit of low saturation our 
simulations yield the expected number of atoms in all four cases. 
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Fig. 3. Simulated single-shiot absorption images for optimum exposure parameters, includ- 
ing fringe removal. The four subplots correspond to the four cases described in section[3] 
Bright colors correspond to high optical density. Each subplot shows 5x5 pixels, with the 
atom initially located at the center. At larger distances from the initial position the PSF 
becomes negligibly small. Using Eq. |9]we extract N=(1.32, 1.39, 0.56, 1.26) atoms for 
subplots a)-d) respectively. These atom numbers are normalized by A'app for each imaging 
setup. 

We find optimum imaging parameters by numerically maximizing the signal-to-noise ratio 
of the central pixel as a function of exposure time and intensity, i.e. maximizing SNR(O) where 

SNR(x) =nW/<7„w (10) 

using the precalculated trajectories of an atom in the integrations of Eqs. (|2]). To avoid con- 
fusion in the case of reflective imaging where the intensity is spatially varying, the intensity 
is here always given in units of saturation intensity for the incoming beam (cf. Fig. [TJ. The 
local intensity at the initial position of the atoms is therefore higher by a factor of four for the 
standing wave. We also determine the SNR of an estimator achieving the Cramer-Rao bound by 
weighting the signal by the mode function of the atomic distribution as described above. SNR 
as a function of exposure time and intensity is depicted in Fig. [2] Table[T|lists the signal-to noise 
ratio for optimum imaging parameters. Note that this can be further improved by a factor « -\/2 
by using the fringe- removal algorithm [18] which reduces the shot noise in the light field by 
optimal averaging of reference images even in the absence of any fringes. Fig.[3]finally shows 
simulated single-shot absorption images of a single atom for the parameters given above, in- 
cluding the use of fringe removal. The number of atoms extracted from these images using the 
optimum estimator is given in the figure caption. 

Depending on the imaging method the effect of recoil blurring on the optical resolution of the 
imaging system can be significant. For our optimum exposure parameters we find the strongest 
recoil blurring in case c), where the RMS position of the atoms in the object plane increases by 
0.9 |xm during the imaging pulse. For all other cases considered here recoil blurring is smaller 
by at least a factor of 3. This bluiiing should be compared to the base resolution of the imaging 
system of 1.2 [im (Rayleigh criterion). Furthermore the images are discretized by the effective 
pixel size of 1 .3 \xm in object space. 



4. Discussion & Conclusion 

In this paper we have outlined the theoretical prospects of detecting single rubidium atoms us- 
ing absorption imaging. We find the optical density of a single atom in single-pass imaging to 
be ultimately determined only by the numerical aperture of the imaging system with a peak 
optical density of OD — 3/2NA^ (l +3/8NA^) in the low saturation limit. We have calculated 
optimum imaging parameters for a number of different cases and could show that the use of a 
reflective imaging setup can significantly improve the signal-to-noise ratio of absorption imag- 
ing. We expect an accuracy of about 82% for the detection of single atoms for the optimum 
parameters given above in case d), which can be improved to 94% (with a false-positive rate 
due to zero atoms of 4%) using the reduction in shot noise achieved by the fringe removal algo- 
rithm 120 1 . This compares to a detection accuracy of only 74% (with a false-positive rate due to 
zero atoms of 15%) for single-pass imaging with the same imaging system. For measuring sta- 
tistical averages in a regular array of traps further improvements can be made using correlation 
analysis |20|. 

We would like to thank Jasper Reinders for help with the simulations and N.J. van Druten 
for helpful discussions of the manuscript. This work is part of the research programme of the 
Foundation for Fundamental Research on Matter (FOM), which is part of the Netherlands Or- 
ganisation for Scientific Research (NWO). 

A. Absorption Signal of a Single Atom 

We start from the simple configuration of a single atom at rest, placed in the centre of the 
minimum waist of a travelling-wave laser beam, see Fig.fTTa). The laser is taken as monochro- 
matte, E.„(.„,,„,.,0 = i4°)(.„,,„)u.(-™) +c.c., where c.c. denotes complex conjugate, u 

is a complex polarization unit vector, (;(:„, y„) is the amplitude in the object plane with 
coordinates x„, yo- We assume that the laser is coUimated in the object plane, so we can 
choose <S'l^\xo,yo) to be real. The incident intensity is ll°\xo,yo) = [<S'ln\xo,yo)]^ /2Zq, with 
Zo = (eoc)"' = ^Mo/eo- 

We describe the total light field as a sum of this incident plane wave plus a scattered wave. 
Assuming they have the same polarisation we ignore from now on the polarization vector The 
total light field is imaged onto a detector plane where we measure the total intensity /det(p = 
(jc,y)). The scattered wave has both a coherent and an incoherent contribution, which we add 
separately as follows, 

/det(p) = ;^|^in(p) + 4r''(p)P+/ir'*'(p) (11) 

/Zo 

= /m(p)+/sc(p)+/if(p) (12) 

The total scattered intensity /sc(p) is the sum of a coherent and an incoherent contribution, 
^^(mcoh)^p^ = slsc°^\p) where s is the saturation parameter ETI . The third term is the interfer- 
ence term, 

/if(p) = ^Re(4l^°''(p)) (13) 

where S'j^°^^ (p) is complex-valued. The interference intensity is typically negative at the center 
of the image, as it reduces the total intensity in the image plane. 
The absorbed intensity is given by 



4bs(p) =/in(p) -/det(p) = -/if(p) -/sc(p) 



(14) 



The terms on the right can be related to the total scattered power fjc after integration. The imag- 
ing system collects a fraction % of the scattered power, while a fraction 1 — is not collected. 
The latter represents the power that is "missing" in the detection plane, by virtue of energy 
conservation: 

Zne = / /sc(p)^/^P (15) 
(l-Z)^sc = / [/in(p)-/det(p)] d^P (16) 



where the integration is over the detection plane. These expressions presume that the incident 
laser beam is entirely captured by the imaging system, so that no interference occurs in the 
noncollected light. When we add the scattered fractions only the interference integral remains. 



Psc = - J Iifip )d'p (17) 

The interference integral must clearly be negative, and equal to minus the total scattered power 
(coherent plus incoherent). Note that it is independent of the collection solid angle. 

Next, we introduce the point spread function to describe how the point-like atom is imaged 
in the detection plane: 

4"''(P)=^/^(P) 

with complex amplitude £/. We choose the (complex valued) point spread function (PSF) to be 
normalized as / p{p)d^p = 1 so that p{p) has the dimension 1/area and a = [J \p{p)\'^ d'^p] ' 
can be interpreted as the effective area over which the scattered light is distributed in the detec- 
tion plane. 

The absorption signal, Eq. ( [14| |, can now be written as 

4bs(p) = Kp(p) - ^ |p(p)|2 (18) 

In most practical cases we can approximate the probe beam as a nearly plane wave, and p{p) as 
sharply peaked compared to (^in(p). Taking the latter outside the interference integral, Eq. ( fTT] ), 
replacing S'i„{p) « (Sin(O), we obtain, using Eq. ( [T3| ), 

//„(P)<<^P^^=»^^-/.. (.9, 
Using the integrals defined above, we can then express the interference term as 



hfip) 



(20) 



If the atom is in the focal plane of the imaging system we can take p{p) to be real-valued 
so that the lm{p{p)) term vanishes. The term also vanishes if the probe beam is exactly on 
resonance with the atomic transition. In this case Im(/9(p)) /Re {p{p)) = tan(argj2/) = 0. The 
approximation of a resonant probe beam is valid as all Doppler shifts are less than 2% of the 
natural Hnewidth of the transition for the optimum exposure parameters found in the main part 
of the paper. In both cases the absorption signal then reduces to 



hUp)~P^c [^t{p{p))-ax\p{p)?] 



(21) 



The leading term in the absorption is simply given by the scattered power. The correction ex- 
presses the partial "refilling" of the absorption dip by scattered light. Within the above approx- 
imation the absorption signal can also be expressed in terms of the scattering cross section, 
using Psc = (7/in(0)/(l +s), 

^ = ^[Re«p))-„|p(p)P] (22, 



